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We study nonequilibrium transport and noise in a generic dissipative 
tight-binding model. Within a real-time path integral approach, we derive 
formally exact series expressions in the number of tunneling events for the 
noise valid for arbitrary bias, frequency and temperature. At zero tempera- 
ture, the low-frequency noise can be summed in analytic form. The resulting 
Shiba-like \u\ -singularity is a consequence of the l/t 2 decay law for the cur- 
rent correlation function. At finite temperature, this singularity is smoothed 
out. 



Low-energy excitations of ID conduction electrons and of edge currents in fractional 
quantum Hall (FQH) systems [IJ are bosonic and constitute prime realizations of a Luttinger 
liquid (LL) phase of interacting fermions [@,|J. A sensitive probe of this collective state is the 
tunneling conductance of electrons through an impurity or point contact in a ID quantum 
wire [|J , and of Laughlin quasiparticles through a narrow constriction for FQH edge currents 
||. In these systems, recent attention has been paid to current fluctuations with emphasis on 
dc noise. To second order in the tunneling amplitude, the dc noise S at T = is classical shot 
noise, S = 2e*I Thus upon measuring dc current and dc noise in the weak-tunneling 

limit, one can determine the charge e* of the tunneling entity. In this way, the fractional 
charge of the Laughlin quasiparticles has been manifested |8|||. In the strong tunneling 
regime, higher order tunneling processes cause a deviation of the dc noise from shot noise 
njll]. Furthermore, correlations between tunneling events lead to colored noise. Singular 



behavior S(u) oc \u\ has been found within a nonperturbative analysis of the equivalent 
integrable Boundary Sine-Gordon (BSG) model [Jl2] , |l3|| . In this treatise, the prefactor of 
the |a>|-law depends on the particular treatment of the reflection matrix for quasiparticle 
scattering at the boundary coming from the Fermi sea (cf. the debate in Refs. |I2| , |I5| ). 

In this paper, we study the low-frequency noise upon employing real-time path integral 
techniques to a generic quantum transport model introduced by Schmid ||T4|| , 



H = H TB + 2J 
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The Hamiltonian (|l|) describes a particle which moves in the tilted tight-binding (TB) lat- 
tice Htb = — f^A Z~2 n ( c n c n+i + c^ +1 c n ) — heq/qo with lattice spacing q , tunneling matrix 
element A, and biasing force F = he/q . The particle is bilinearly coupled via its position 
Q = Qo J2 n c li c n to a bath of delocalized harmonic bosons. In the reduced density matrix 
description for the Brownian quantum particle, all effects of the environmental coupling are 
captured by the spectral density J(u) = l^ a {c 2 a /m a uj'^)5{uj — uj a ). The TB model 

(0) is dual to a corresponding weak-binding (WB) model, in which the particle moves in 
a tilted cosine potential. In the duality, the power s in the generic form J(cj) oc uj s is 
transformed into 2 — s. Thus, super-Ohmic friction (2 > s > 1) is mapped on sub-Ohmic 



friction (0 < s < 1), and vice versa [15]. In the Ohmic case, J{ui) = 2Kuj, the TB and WB 
models are selfdual to each other. The dimensionless Kondo parameter K is mapped onto 
the reciprocal value 1/K, and the bias e onto e/K, as first observed by Schmid [[14]. 



The TB and WB model describe many transport problems in solid state physics. Re- 
garding applications of the model ([!]) to the beforementioned case of tunneling of a charge 
in a LL phase, the spectral density is strictly Ohmic for u ^ u c . In the scaling limit, uj c is by 
far the largest frequency of the problem. The interaction parameter g of the Luttinger model 
corresponds to 1/K, and the bias in the TB model (|l|) is related to the voltage drop V across 
the weak link or constriction by he = eV, where e is the unit charge. Another important 
example where the model (|Ij) applies is a voltage-biased Josephson junction [|1^,|T^]. Here, 
the WB and TB limits correspond to the charge and phase representation, respectively. 

Fluctuations are measured relative to a stationary state. Therefore, we proceed as fol- 
lows. We first presuppose that the global system is in a factorized system-bath state in the 
infinite past. Then we switch on the bath coupling and the tilting force, and let the system 
evolve under reign of the full Hamiltonian ([!]). By the time t — 0, the system has reached 
a stationary nonequilibrium state, indicated by the index "ne". The fluctuations of the 
current j(t) = eq(t)/qo are reflected in the noise spectrum, which is the Fourier transform 
of the symmetrized current correlation function, 

/oo 
du;e^(j(t)j(0)+j(0)j(t)) ne . (2) 
-oo 

Observing that the current correlation function is related to the mean square displacement 
C nc (t) = ([q(t)-q(0)] 2 ) nc hy 

(j(t)j(0)+j(0)j(t)) ne = e ^C nc (t) , (3) 

we can express the noise S(u>) in terms of the spectral function of C ne (t) as 

e 2 ~ ~ f°° 

S(u) = --u 2 C nc (u) , C nc (u)= dte^C ne (t) . (4) 

% J -co 

In particular parameter regimes, the noise can directly be expressed in terms of transport 
quantities, e. g. the shot noise S = 2eVG(V) and Johnson- Nyquist noise S = ik^TGiV = 
0), where G(V) = I(V)/V is the nonlinear conductance. For later convenience, we also 
introduce the differential conductance G&(V) = dI(V)/dV. The conductances are related 
to the mobilities of the model (P, fi(F) = (q(t — > oo)) fc /F and f^d(F) = d(q(t — > oa))f c /dF, 
by G = e 2 /i/gg and G^ = e^n&jq 1 . The index "fc" indicates that we may choose for the 
system-bath state at initial time zero a factorized form. This is possible since the leading 
contribution to mean values at asymptotic times t — » oo does not depend on system-bath 
correlations in the initial state. 

Within nonequilibrium linear response theory, the nonequilibrium susceptibility 

Xne(t) = T<3>(t) (q(t)q(0) - q(0)q(t)) ne , XnoH = / dte^ Xne (t) (5) 
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comprises the differential mobility as 

Hd{F) = lim [ - iuxnci^)] = lim Xne(t) • (6) 

u;^0 t^oo 

Formally exact series expressions in A for the above quantities can be derived upon 
employing a real-time path integral method for the reduced density matrix (RDM). Because 
of the discreteness of the TB model, a simple picture of a grandcanonical gas of charges 
can be developed. The charges are stringed at the flip times {tj}, and the effects of the 
environment are in the charge interactions. For a path with n flips in the negative-time 
branch (ti, t 2 , ■ ■ ■ , t n < 0) and m flips in the positive-time branch (t n+ i, t n+2 , . . . , t n+m > 0) 
labelled by charges = ±1}, {rjj = ±1}, the interaction factor or influence function is 

Fn,m Gn,mHn,mi where 



■ n+m k — 1 s. f n+m—1 ^ 

^2 ^'( tfe _ *i) f ' Hn > m = exp y Vflj^+m > (7) 

• k=2 j=l ' ^ j=l ' 



with 6j^ n+m = Y^k=j+\ £kQ"(tk — tj). The charges {r]j} describe propagation on the infinite 
RDM in diagonal direction, whereas the charges indicate moves perpendicular to the 
diagonal, i.e. represent quantum fluctuations. The pair interaction Q(t) is the second 
integral of the autocorrelation function of the random force f(t) = ^2 a c a x a (t) in thermal 
equilibrium, Q(t) = (qo/h) 2 (f(t)f(0)) eci In the Ohmic scaling limit, the pair interaction 
reads 

Q(t) = Q'(t) + iQ"{t) = 2K In [(h/3u e /ir) sinh(vr|t|/^)] + iirK sgn(t) . (8) 
The influence of the bias e or voltage drop V is included in the bias phase factor 

n+m n+m— I 

K,m = exp(-20„ jm ) , </y m = e 22 &i = e ^2 9l,n+mPl ■ (9) 

1=1 1=1 

Here, we have introduced as an off-diagonal measure the cumulative charge after I jumps 

gi, n +m = Y2=l+i €k = - Y,k=i &j an d tne time intervals p t = t w - U. 

Following the lines of Ref. flI5 |, it is straightforward to arrive at formally exact series 
expressions for the Laplace transforms of the above quantities in which the summation over 
all possible sequences {r/j} has been performed already. 

The first moment in the stationary nonequilibrium state with initial condition (q(t = 
0)) ne = is given by 

(qW)ne = (<?(A)) fc + i? ne (A) ■ (10) 

The term (g(A))f c is the first moment with factorized system-bath initial state at time zero. 
This term has dynamics only in the positive-time branch. The residual term i? ne (A) describes 
the corrections due to the correlated initial state, 

^ 00 POO 

(9(A))fc = ^2 ^2 / ^ , o,m(A)^ao^G o , m sin0 Oim , (11) 

mcvcn^ 

^ 00 POO POO 

R ne (X) = — ^ J V n ,M) J dl ~n ds n e' Xs " ^2 a nin G n,m SHI n , m . (12) 



n,m=l'" J ,/u {&} 



3 



The sum is over all charge sequences which form a neutral cluster, XwSl™^ = 0- The 
integration symbol includes the A-factors associated with each jump and the integrations 
over the interval lengths except for the interval p n , 



poo / n + m 1 poo \ n 1 i*oo 

/ Pn,m(A) ••• = A n+m ( J] / d Pz e- A " HZ ( 13 ) 



Z=ri+1 uyJ / k=l 



The interval p n around t = is divided into the negative-time part r n = —t n and the 
positive-time part s n = t n+ i and is treated separately. Finally, the phase factor 



n+m—l 

„,„-'.,i l) in+m - 2)/2 J] «in(^ >+m ) (14) 



subsumes the interaction between the {rjj} and charges. In the Ohmic scaling limit, 
the phase 0j t n+ m is independent of the flip times and is 6j^ n+m = nKg^ n + m . 
Similarly, the series expression for the nonequilibrium susceptibility reads 



oo 

Xnc(A) = TTT 

iriX 

n, 



°° POO POO 

^ / £>n, m (A) / dr w ds n e~ As "y^ g n , n + m a$ m G n ,m cos n ,, 



The differential mobility can be obtained either from this form using the relation (Q) or from 
its definition with eq. (pTTj) . In any event, we arrive at the series expression 



2 oo „oo^ sm-i \ 

^=11 J2^ m/2 ~ 1 ^ a -> ( n sin °i,m)G , m B 0)m , (is) 

m even ^ {^} ^ Z=l ' 

^o,m = ^- sin o , m = ( ^z, m p/ J cos o , m • 

As well, the mean square displacement can be divided into the second moment with 
factorized initial state at time zero, (g 2 (A))f c , and the correction term K ne (X), 

(7 nc (A) = (g 2 (A)) fc + it nc (A), (16) 

00 POO 

(<f (A))fc = ^ ^ / ^0,m(A) ^4^Go,mCOS0 O ,m , (17) 
m even ® { } 

^ne(A) = - ^ / T> n>m (A) dr n ds n e~ Xsn ^ a^lnG^rn cos <p ntm . (18) 

A n,m=l J JO te} 

The phase factor resulting from the {r/j} summation differs from ctn,L- We have 

n+m—l 

«Sn = -iQO a£m XI Cot (^,n+m) • (19) 
fc=n+l 
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The formally exact series expression for the noise follows from (f|) and ([H]) - fll8|) as 

S{u>) = 2(e 2 /g 2 ) ^ [A 2 (7 nc (A)] A= _ +0+ . (20) 
Consider now the A-expansions of eqs. (|10D - (|i~2D and eqs. ( fTED - fllSf) about A = 0, which 



determine upon Laplace transformation the behavior of the respective quantities at long 
times. We first observe that the phase factor cinjn vanishes if the system reaches a diagonal 
state gi^n+m = 0. Therefore, all paths that contribute to the first moment hit the diagonal 
of the RDM only in the initial and final state. Thus each interval pi separates two charged 
clusters, and for A = 0, the interaction factor G ntTn at T > or the bias factor B nm at T = 
ensures the convergence of the respective integration. For the first moment with factorized 
initial conditions, eq. (|TT|), we find 

A 2 (g(A)) fc = F/i + \ qoo + 0(X 2 ) , goo = \im(d/d\) [A 2 (g(A)) fc ] . (21) 

The series expression for p is as in (|15D , but with B 0tm = sin(0 o ,m)/e- 

The series (0) can be rearranged upon using the relation (we put p n = r n + s n ) 

/ dr n ds n e- Xsn A{r n + s n ) = - / dp n [ 1 - e~ xpn ] A(p n ) , (22) 
Jo A Jo 

where A(p n ) subsumes all correlations across the intervall p n . It is then straightforward to 
see that the A-dependence of the expression A 2 i? ne (A) is fully cancelled by A 2 (g(A))f c . In the 
end, we find using eq. (|10D for all A 

\ 2 (q(\)) nc = Fp. (23) 

This verifies that the mean position with initial value (q(t = 0)} ne = evolves in the 
stationary nonequilibrium state for all times as (q(t)) ne = Fpt. 

After these preliminaries, consider now the mean square displacement, eq. (p~6|) . In 

(2) 

the phase factor On,m, the factor cos(7r Kg^ n+rn ) can be allocated to any interval in the 
positive-time branch. However, upon employing the relation fl22f) , the series for (7 ne (A) = 
(g 2 (A))f c + i^ nc (A) can be rearranged so that most contributions of the sum in eq. (0) cancel 
out. The only remaining contribution is the one in which the cosine factor covers the first 
interval in the positive-time branch, labelled by p n . Thus, this is the only interval which 
is allowed to be occupied by a diagonal state, g n)n+m = 0. Since p n separates two neutral 
{£j}-clusters, the attractive bias and charge interaction factors are missing, and hence this 
interval gets very long as A — > 0. It is convenient to factorize the interaction term G n ^ m (n, m 
even) into the intra-cluster interaction factors G n fl and Go,m and the inter-cluster interaction 
factor G in t, G n , m = G n flGo^ m G m t. Expanding the inter-cluster interaction stretching over 
the long interval p n , we get 

n— 1 n+m— 1 

9l,n+mPl / j 9k,n+mPk ■ 

(24) 

1=1 k=n+l 

The first term is the contribution in the absence of the cluster interaction, and the second 
term represents the dipole-dipole interaction. Accordingly, we get two contributions to 
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C ne (X), called die 1 (A) and C®(\). In the first case, the cluster to the left of p n is just 
Fp = A 2 (g(A)) ne , whereas the cluster to the right of p n represents A 2 (g(A))f c [Note that each 
cluster is actually an infinite series in A 2 ]. Since the two clusters are noninteracting, the 
integration over p n gives a factor 1/A. Thus we find 

A 2 (7^(A) = 2[Fp] [A 2 (g(A)) fc ]/A = 2F 2 p 2 /\ + 2F Moo + 0(A) . (25) 

In the second form, we have used the result ([H]). 

Consider next the dipole contribution Cj$(\). Upon inserting the second term of eq. (p4|) 
into the series for C ne (X), we observe that the two sums can be generated by suitable differ- 
entiations with respect to the bias. Thus we find 

2 00 00 poo poo 

A 2 C£HA) = -^E E / dp n e-^Q'(p n ) V n , m (X) (26) 

n even m even ^ " 

^2 E G! n,o G! o,m sin <f> nJ aj (— sin </> , m J JJ sin0j >n+m . 



The sums S{&}„ an d S{£ fe } m cover & U neutral charge clusters, i.e., they are constrained by 

Er =1 6 = oandErr + i6 fc =o. 

At T = 0, the dipole-dipole-interaction is found from eq. (|j) as Q'(t) = —2K/t 2 . Sub- 
stituting this form, the integration over p n in the limit A — > yields 

POO 

/ dp n e~ x ^ Q'(p n ) « — 2KT A ln(A/A*) . (27) 

Now, in leading order of A both the cluster on the left and the cluster on the right of p n can 
be identified with the differential mobility pd, eq. (|l5|) . Thus we find 

A 2 ^ (A) = - (^ 2 K/q 2 )p 2 d A ln(A/A*) [1 + O(X)] . (28) 

So far, we have considered the contributions to C nc (X) in which the period p n the system 
dwells in a diagonal state is very long. However, there are also contributions in which the 
diagonal state is either occupied for a short period only or not occupied at all, g n ,n+m 7^ 0. 
In both cases, the charges form a single neutral cluster with regular behavior in the limit 
A — > 0. These terms are combined in the expression 

X 2 C^(X) = 2D' + (D(X). (29) 

Putting all contributions together, the mean square displacement takes the form 

2 Ati 2 K 

A 2 C nc (A) = —F 2 p 2 + (2Fp goo + 2D') —p\ A ln(A/A*) + 0(A) . (30) 

A % 

Turning to the time regime, the expression ( |30"D gives the behavior at long times, 

C ne (t) = F 2 p 2 t 2 + 2Dt + (4k 2 K/q 2 ) p\ \n{t/t*) . (31) 
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Here, all terms linear in t are combined in the diffusion constant 



D = ~ hm 1 [C nc (t) - (q(t))l] =D' + F Woo . (32) 
The current correlation function at long times is found from eqs. (|3|) and (|3~ID as 

(i(0i(0) +j(0)j(t)U = 2(eF^/ qo y - M . (33) 



Here, we have normalized the differential mobility on the mobility of a free Brownian particle, 
fi& = /-td//-to with A*o — qo/2iihK. 

Using (Effi), it is easy to commute the A-expansion ( |30|) into the low-frequency noise 
spectrum. The 1 / A-term in the expansion (|30|) does not contribute to the noise. The term 
of order A directly relates the dc noise S(u = 0) to the diffusion constant, 

S(uj = 0) = 4e 2 D/q 2 . (34) 

In the integrable BSG model, a basis of interacting quasiparticles which scatter off the 
impurity or point contact one by one can be introduced. Upon generalizing a Landauer-type 
approach familiar from transport of free electrons to this interacting case, it has been argued 
that the dc noise at T = is connected with the conductance by differentiation with respect 
to the voltage [10| or the tunneling amplitude []Tl,20| 



S(u = 0) = eV A dG/dA = [1/(K - l)]eV 2 dG/dV . (35) 

In order A 2 , this formula reproduces shot noise, S(u = 0) = 2eVG = 2el. The higher order 
terms describe corrections due to quantum coherence for strong tunneling. The relation ([35]) 
can be proven on the basis of (|TJ) - ( pT8|) and is a consequence of (i) a generally valid detailed 
balance relation (DBR) that follows from analytic properties of the pair interaction Q(t) and 
(ii) of a special DBR which is restricted to the scaling limit, i.e. relies on the particular 
form (|8]) for Q(t) [[Hp. The general DBR states that at T = there is no real occupation of 
TB states higher than the initial state. The special DBR states that only those paths for 
which all {rjj} charges have equal sign, i.e., which are free of backward hops parallel to the 
diagonal, contribute to the noise. 

The low-frequency correction AS(uj) is captured by the term of order Aln(A/A*) in 
Using (|20|), we obtain for the noise the nonperturbative universal form 



ASM = 2^ (e7<Zo) Jr d M = 2hG (g)G 2 d \u\ . (36) 

In the second form, we have used transport quantities of the related Luttinger model, i. e. 
the dimensionless differential conductance G*d — Gd/Go(g), normalized with the static limit 
of the microwave conductance of a quantum wire or the Hall conductance, Go(g) = ge 2 /h. 



This result agrees with the findings of Ref. [13 



The algebraic 1/t decay law of the current correlation function (|33|) and the resulting 
| uj | -singularity in the noise spectrum ( |3T)D are signatures of the unscreened dipole-dipole 
interaction. These forms are universal since they do not depend explicitly on the parameter 
K or g and the applied bias. The parameters are only in the mobilities or conductances which 
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determine the prefactor of the 1/t 2 - and [a; [-law. The forms (^) and fl36|) are analogous to 
the 1/t 2 decay law of the position correlation function [18|,19| in the Ohmic two-state system 



and the equivalent Shiba relation in frequency space, respectively |18|,|T9[. The difference, 
however, is in the prefactor. For the position autocorrelation function, each of the two sets 
of dipoles corresponds to the static susceptibility xo, whereas in the present case each dipole 
set represents the differential mobility, which is related to the nonequilibrium susceptibility 
by eq. (^). In fact, in transport problems with a stationary nonequilibrium state, the static 
susceptibility diverges, and the differential mobility (id{F) is the proper quantity here. 
At finite temperature, the dipole-dipole interaction is given by 

Q\t) = -2K {7c/h(3) 2 / sinh 2 (vrt/^) . (37) 

Performing the integration over the interval p n between the dipoles, eq. (0), and switching 
to the noise, the low-frequency contribution from A 2 Cl^(A) takes the form 

A5( Wl T) = [Aj 2 u coth (h^/2) « [fi d ] 2 + ^ [fc] 2 u 2 + 0{^) . 

Thus, the | a; [-singularity is smoothed out at finite temperature. We get a contribution to 
the dc noise which is proportional to T and a correction term of order uj 2 . However, there 
are additional contributions oc u and oc lu 2 from \ 2 C£\\) and \ 2 C$(\). However, the 
corresponding prefactors can not be given in analytic form for all K. 
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